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Abstract 

We investigate the ground-state property of a one-dimensional two-orbital Hubbard model at quarter filling by numerical techniques 
such as the density-matrix renormalization group method and the exact diagonalization. When the Hund's rule coupling J is zero, 
the model is SU(4) symmetric. In fact, both spin and orbital correlations have a peak at q=ir/2, indicating an SU(4) singlet state 
with a four-site periodicity. On the other hand, with increasing J, it is found that the peak position of the orbital correlation 
changes to q=~n, while that of the spin correlation remains at q—n/2. We briefly discuss how the SU(4) symmetry is broken by J. 
© 2008 Elsevier B.V. All rights reserved. 
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It has been widely recognized that the interplay of spin 
and orbital degrees of freedom plays a significant role in the 
emergence of exotic magnetism in strongly correlated elec- 
tron systems with orbital degeneracy. In this context, the 
highest symmetric SU(4) spin-orbital model has been one 
of the subjects of much interests from a theoretical view- 
point. In particular, the one-dimensional model is Bcthc 
ansatz solvable [1], and the combined quantum effects of 
spin and orbital have been revealed by analytical and nu- 
merical investigations [2,3,4]. Indeed, it is characteristic of 
an SU(4) singlet state that correlation functions show criti- 
cal behavior with a four-site periodicity and the elementaly 
excitation is gapless. 

The highest SU(4) symmetry originates in the situation, 
where electrons hop only between the same types of orbitals 
with equal amplitude and the Hund's rule coupling is ig- 
nored in a two-orbital Hubbard model. In a more realistic 
situation, however, the Hund's rule coupling should break 
the SU(4) symmetry down to SU(2) spin xU(l) or bitai [3,5,6]. 
To clarify the effect of such symmetry breaking, Lee et al. 
have studied an SU(4) Hubbard model perturbed by the 
Hund's rule coupling by means of renormalization-group 
and bozonization methods, and proposed that the spin gap 
opens for an arbitrarily small Hund's rule coupling [5]. 

In this paper, we investigate spin and orbital correlations 
in a one-dimensional two-orbital Hubbard model with one 



electron per site. The Hamiltonian is given by 

H = (4ra d i+ira + h.C.) + U ^ Pir]Pirl 

i,T.<j i,r 
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where di T<T is the annihilation operator for an electron with 
spin a in orbital r(=a, (3) at site i, pi ra =d\ ra di rcr ^ and t is 
the hopping amplitude. Hereafter, t is taken as the energy 
unit. U , U', J, and J' denote intra-orbital Coulomb, inter- 
orbital Coulomb, exchange (Hund's rule coupling), and pair 
hopping interactions, respectively. Note that U =U'+J+J' 
due to the rotational invariance in the local orbital space, 
and J=J' due to the reality of the wavefunction [7]. When 
the Hund's rule coupling is zero, i.e., U~U' and J=J'=0, 
the system possesses the SU(4) symmetry. 

We analyze the model with 32 sites (7V=32) in the open 
boundary condition by the density-matrix renormalization 
group method [8] . The finite-system algorithm is employed 
with keeping 400 states per block and the truncation error 
is estimated to be 5 x 10 -6 at most. We also use the Lanczos 
method for the analysis of a four-site periodic chain. In this 
paper, we set U'—J=20 to consider the strong-coupling 
region, and investigate the dependence on J. 

In Fig. 1, we show our DMRG results for the spin and 
orbital structure factors, defined by 
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Fig. 1. Spin and orbital structure factors for (a) J=0, (b) J=0.5, 
and (c) J=l. 

S(q) = Y,j q(J - k) {S!S£)/N, (2) 

T^^^^HT^/N, ( 3 ) 

with S , f=X)r(/ ^T-PiTi)/2 &ndTf=J2 a (piaa-Pi0cr)/2- At 
J=0, 5(g) and T(q) coincide with each other, and we can 
observe a peak at q—Tr/2, clearly indicating an SU(4) sin- 
glet state with a four-site periodicity, which is consistent 
with the previous numerical work for the SU(4) spin-orbital 
model [3]. On the other hand, 5(g) and T(g) exhibit dis- 
tinct behavior with increasing J, as shown in Figs. 1(b) and 
1(c). It is observed that the peak position of T(q) changes 
from q=n/2 to g=7r, since T(w/2) is suppressed, while T(ir) 
is enhanced due to the effect of J. As for the spin correla- 
tion, the peak position of 5(g) remains at g=7r/2 even for 
finite values of J, since the magnitude of S(n/2) increases 
in sharp contrast to the case of T(q). 

To obtain intuitive understanding for the changes in 5(g) 
and T(g), it is useful to consider a four-site system, which 
is a minimal model to form the SU(4) singlet ground state 
at J=0 [9]. The SU(4) singlet state is expressed as 

|SU(4)> = (l/x/24) Y, 4 T 4l d U d wl0), (4) 

where |0) is the vacuum state and the summation is taken 
over all permutations for site indices. Note that |SU(4)) 
consists of 24 states with the same weight for each state. 
For finite J, however, the ground state is not represented by 
|SU(4)) itself. As shown in Figs. 2(a)-(c), the 24 states are 
split into three classes with eight states for each according 
to the weight in the spin-singlet ground state [6] . Note that 
each class is characterized by the peak positions of 5(g) 
and T(g), denoted by (g sp i n , gorbitai): (7r/2, 7r) for the class 
a, (ir/2, 7r/2) for the class b, and (n, ir/2) for the class c. 

In Fig. 2(d), we show the J dependence of the weight of 
each class m in the ground state, 

w m = Y\(<f>i\i>G)\ 2 , (5) 

where Vg is the ground state and fa denotes the basis. 
At J=0, three classes contribute to the ground state with 
equal weight. With increasing J, the weight of the class a 
increases, while those of the class b and c decrease and the 



Fig. 2. (a)-(c) Three classes of electron configuration in the SU(4) 
singlet state of the four-site system. Each class has eight equivalent 
states due to the translation and the change of orbitals. (d) Weight 
in the ground state of each class in Figs. 2(a)-(c). 

total weight of the three classes does not change. In the class 
a, as shown in Fig. 2(a), partly spin ferromagnetic (FM) 
alignment appears due to the Hund's rule coupling, and an 
antifcrro-orbital (AFO) configuration is favored to avoid 
the energy loss in the hopping process due to the intra- 
orbital Coulomb interaction, indicating the instability to a 
FM/AFO state due to the Hund's rule coupling. Thus, the 
correlations of S(n/2) and T(n) are enhanced, leading to 
the change of the peak position of T(q), as shown in Fig. 1. 
We note that with further increasing J, the ground state is 
changed to a FM/AFO state. 

We have tried to understand how the spin gap A s of 
the present model (1) is affected by J, but it is a difficult 
task to estimate A s with precision, because of charge, spin, 
and orbital degrees of freedom. To clarify the behavior of 
A s in the thermodynamic limit, it would be appropriate 
to investigate an effective spin-orbital model in the strong- 
coupling limit, which is an interesting future issue. 

In summary, at J=0, 5(g) and T(q) agree with each other 
and have a peak at g=7r/2 due to the SU(4) singlet ground 
state. On the other hand, with increasing J, we observe the 
transition of the peak position of T(q) to g=7r in accordance 
with the change of relevant spin-orbital configuration. 

We thank K. Ueda for useful discussions. T.H. is sup- 
ported by the Japan Society for the Promotion of Science 
and by the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan. 
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